Wave packet dynamics in 2DEG with spin orbit coupling: splitting and zitterbewegung 
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We study the effect of splitting and zitterbewegung of ID and 2D electron wave packets in the 
semiconductor quantum well under the influence of the Rashba spin orbit coupling. Results of 
our investigations show that the spin orbit interaction induces dramatic qualitative changes in the 
evolution of spin polarized wave packet. The initial wave packet splits into two parts with different 
spin polarization propagating with unequal group velocity. This splitting appears due to the presence 
of two branches of electron spectrum corresponding to the stationary states with different chirality. 
It is demonstrated also that in the presence of external magnetic field B perpendicular to the electron 
gas plane the wave packet splits into two parts which rotates with different cyclotron frequencies. 
It was shown that after some periods the electron density distributes around cyclotron orbit and 
the motion acquire an irregular character. Our calculations were made for both cases of weak and 
strong spin orbit coupling. 

PACS numbers: 73.21. Hb, 71.10.Pm, 72.10.-d, 73.23.-b 
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I. INTRODUCTION 

Producing and detecting spin polarized currents in 
semiconductor nonmagnetic devices is the ultimate goal 
of spintronics. The intrinsic spin orbit interactiorji ex- 
isting in low dimensional systems which couples electron 
momentum to its spin is one of the most promising tools 
for realizing spin polarized transport. For this reasons, 
during the last years a substantial amount of work has 
been devoted to study effects of spin orbit interaction on 
the transport properties of nanostructures (for a review, 
see, e. gj^i^ii). 

At first time the electron wave packet dynamics includ- 
ing the problem of zitterbewegung in semiconductor quan- 
tum well under the influence of the spin orbit Rashba and 
Dresselhaus coupling has been considered by Schliemann, 
Loss and Westervelt^. In this work the oscillatory motion 
of the electron wave packets reminiscent of the zitterbe- 
wegung of relativistic electrons was studied for free elec- 
tron motion i.e., in the absence of electric or magnetic 
field. The authors oS^ predicted the resonance amplifica- 
tion of zitterbewegung oscillations for the electron moving 
in a quantum wire with parabolic confinement potential 
and propose to observe this fundamental phenomena ex- 
perimentally using high resolution scanning probe mi- 
croscopy imaging techniques. 

The zitterbewegung of the heavy and light holes in 3D 
semiconductors was investigated irt^.Tn this paper the 
semiclassical motion of holes in the presence of a constant 
electric field was studied by numerical solution of the 
Heisenberg equations for momentum and spin operators 
in the Lattinger model of spectrum. It was shown that 
the hole semiclassical trajectories contain rapid small am- 
plitude oscillations reminiscent the zitterbewegung of rel- 
ativistic electrons. It should be noted, however, that the 
spatial structure of the wave packet and the changing of 
its shape due to effect of splitting in^'^ was not consid- 
ered. 



At the same time the splitting of spin polarized elec- 
tron beams in the systems with spin orbit coupling was 
investigated in a series of works. In particular, the au- 
thors of papers^'^ propose to use the lateral interface be- 
tween two regions in gated two-dimensional heterostruc- 
ture with different strength of spin orbit coupling to po- 
larize the electron. They have shown theoretically that in 
this structure a beam with a nonzero angle of incidence 
splits into some spin polarization components propagat- 
ing at different angles. The similar effect of electron spin- 
polarized reflection in heterostructures and spatial sep- 
aration of the electron beams after reflection has been 
observed experimentally in^. 

The transverse electron focusing in systems with spin 
orbit coupling at the presence of perpendicular magnetic 
field was theoretically analyzed ini^ where it was shown 
that in the weak magnetic field regime and for a given en- 
ergy, the two branches of states have different cyclotron 
radii. The effect of spatial separation of the electron tra- 
jectories of different spin states in a perpendicular mag- 
netic field has been experimentally observed in^^. 

In this work we study the striking dynamics of the elec- 
tron wave packets in a narrow ^355 quantum well at the 
presence of the spin orbit fc- linear Rashba coupling, which 
arise due to structural inversion (" up-down" ) asymmetry. 
The splitting of the wave packets in two parts appear due 
to the presence of the electron states with "plus" and 
"minus" chirality, which propagate with different group 
velocity. These two parts of the split packet can be char- 
acterized by different spin density. It is found that elec- 
tron trajectories contain small amplitude damped oscil- 
lation. We show that the packet splitting leads to the 
damping of zitterbewegung. The splitting and zitterbe- 
wegung of wave packet is naturally accompanied by its 
broadening due to effect of dispersion. 

We investigate also the atypical cyclotron dynamics of 
the wave packet in a perpendicular magnetic field. It 
was shown that due to the spin orbit coupling the packet 
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with spin parallel to the magnetic field splits into two 
parts which rotate with different cyclotron frequencies. 
We determine the moments when two parts of the packet 
are located at opposite points of the cyclotron orbit and 
after that they return many times back to their initial 
state. With the time due to the incommensurability of 
the cyclotron frequencies and the ordinary packet broad- 
ening the electron density distributes randomly around 
the cyclotron orbit. All our calculations were made for 
the material parameters of the real semiconductor struc- 
tures with a relatively strong and weak spin orbit and 
Zeeman interaction. 

The paper is organized as follows. In Sec. II we intro- 
duce the Green functions for two dimensional electrons 
in the presence of Rashba spin orbit interaction and ana- 
lyze the evolution of ID wave packet. The analytical and 
numerical results illustrate the effects of packet splitting 
and zitterbewegung. In section III we describe in details 
the time development of the 2D wave packets. Finally, 
in Sec. IV we discus the manifestation of the spin or- 
bit interaction in the evolution of coherent wave packet 
in a magnetic field perpendicular to electron gas plane. 
The splitting of the initial coherent packet and distribu- 
tion electron probability via cyclotron orbit is considered. 
Section V concludes with a discussion of the results. The 
Appendix provides the mathematical details necessary to 
obtain Eqs. (36a) and (36b). 



II. THE DYNAMIC OF THE 
ONE-DIMENSIONAL WAVE PACKETS 

In this section we consider the specific character of the 
wave packet dynamics in the systems with Rashba spin 
orbit coupling^. The Hamiltonian of the system under 
consideration reads 
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(1) 



where p — —iTiV is the momentum operator, m is the 
electron effective mass, a is the Rashba coupling con- 
stant, and the components of the vector a denotes the 
spin Pauli matrices. The eigenfunctions for the in-plane 
motion identified by the quantum numbers 'p{PxTPy) are 
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Here (/? is the angle between the electron momentum p 
and X axis, so e"^ — Ee^^Ejl s — ±1 denotes the branch 
index. The energy spectrum of the Hamiltonian (1) cor- 
responding to two branches has the form 
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To analyze the time evolution of electron the initial 
states we use the Green's function of the nonstationary 
equation, which is a non diagonal 2x2 matrix 
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Here i,k = 1,2 are matrix indexes and matrix elements 
can be written as an integrals 



G,fc(r,r',t) 



dp(/)p,,,(r,i)(/);,,,(r',0). (6) 



In the present section we examine in details the dy- 
namics of the quasi-lD wave packet in 2D system with 
spin orbit coupling. This problem allows the analytical 
solution. Let at the initial time t = wave function to 
be a plane wave with wave number pox modulated by a 
Gaussian profile and spin polarized along the z direction 
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Ly is the size 



where coefficient C is equal to (ji^^) 
of the system in the y direction. The variance of the posi- 
tion operator < (Ax)^ > in this case is equal to ci^/2 and 
the variance < (Ay)^ > exceed this value. The variance 
of the momentum operator px is < (Apx)^ >= /2cP, 
and the average p is equal to pox ■ One may consider the 
initial wave function as the limiting case of a 2D packet 
with the width along y direction much greater than along 
X i.e., Ly ^ d. 

The electron wave function at any arbitrary moment of 
time can be found with the help of the Green's function 



^i{x,t) 
*2(2:,t) 



= / dx'dy' 
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where matrix elements Gn and G21 of the matrix (5) are 
determined by Eqs. (2), (3) and (6) 
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Zm 
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where p — \Jp^+p1- Using the definition v = ^ = 
j^[H, r], one can obtain from Eq.(l) the velocity operator 



G2i(r,r',i) = 



exp(- 



in^t p(r — r')s 
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X sm(— — ) 

n p 
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By using the formula 

n=l 

+2z^J2n-i(<?)sin((2n-l)V'). (11) 
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and by integrating over the angle variable in Eqs.(9), (10) 
we finally have 



Gil = / exp(-i^Jo( ^ )x 



27r/i^ y ' ' 2m;i ' /i 



apt 

X cos[^- )pap, 
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^21 = T. ; / exp(— ^- — -x 

27r?i>-r'| J ^ 2mn 

^ ,p\r- r'l . ,apt. , 
X Ji ( — ) sin ( — )pdp, 
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where Jq and Ji are Bcsscl functions. Substituting Eqs. 
(12), (13) and (7) into Eq.(8) and integrating over x' and 
y' , we find the analytical expression for the spinor com- 
ponents ■)/'i,2(a;, i).It should be noted that two electron 
bands with chirality "plus" and "minus" give different 
contribution to the electron wave functions. The calcu- 
lation of the expressions for ^'1,2 leads to the following 
electron probability densities |^'i|^ and |4'2p at any ar- 
bitrary moment of the time 
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, (x+(a — hko/m)t)^. 

-P( ,2(1+^2,2) ) + 



{x- {a + hko/m)t)'^\ 
^ 2^2(1+72^2) 

,2(ko(f + ^tx)at. 
X cos(^:327— — 5-2^) 
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d2(l+72,2) ) + 
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rf2(l+72t2) )- 

(x + (a - fikn / ■m)t)'^ 
2^2(1 -F72f2) 



(x- {a + fikQ/m)tY\ 
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where 7 = Ti/(Pm is the inverse broadening time pox = 

As follows from Eqs. (14a), (14b) the shape of the 

function p{x^ t) essentially depends on the parameter 

2 2 ,2 

11 = ^2 ■ In the case wide packet when the momen- 
tum variance is much more (ma) 2 and the inequality 
7/ -C 1 takes place, the evolution looks like at the ab- 
sence of Rashba term. Otherwise when 77 1 the initial 
wave packet splits into two parts which propagate with 
different group velocity, so the distance between these 
two parts increases linear in time. This two parts cor- 
respond to the first and second terms in square brackets 
in Eq.(14a) and Eq.(14b).The third terms in Eq.(14a), 
(14b) describe the oscillation of the components of elec- 
tron density |*i|2 and |'I'2p in the region of the overlap- 
ping of two split parts of the packet. It is clear that these 
oscillations originates from the interference between the 
statcis of different spectrum branches. When two parts 
of the packet move away from each over the amplitude 
of the oscillations decreases. The period of these oscilla- 
tions along the x direction depends on the initial width 
of the packet d and equals to Aa; = ^(^{l + ^"^t^) / a'yt^ . 
So, if inequality ^1 takes place the period of oscilla- 
tion decreases with time and equals to Ax = irmd^ / aht^ 
and when ^ 1 the oscillation period is not depend on 
the time Ax — T:1i/ma. 

To illustrate the evolution of the electron probability 
density p(x, y) = |*ip + 1*2|^ we plot this function using 
Eq.(14a), (14b) at the Fig. 1(a) for the moments of the 
time: t = Q,t= 1, 5, i = 7 (in the units of tq = 7~^). The 
calculations was made for the parameters GaAs/InGaAs 
electron system: m = 0, 05mo, a = 3, 6-10^ cm-sec""'^ and 
the packet parameters: d = 10~^ cm, fco = 2, 5-10^ cm~^. 
Here one can clearly see that initial Gaussian wave packet 
Eq. (7) splits up at i > into two parts propagating along 
the X direction. The width of each part of the packet 
increases in time as for the case of free particle. 

To analyze spin dynamics one can consider the time 
evolution of the spin density 



Si{x,y,t) = ^{^l,%)ai 



*i 

*2 
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Using Eqs. (14a) and (14b) we immediately find the ex- 
pression for spin density = |(|*i(r, t)|2 — |\E'2(r, t)|2), 
which demonstrate the oscillatory behavior as a function 
of X (see Fig. 1(b)). The period of oscillation here is the 
same as for the functions |*i,2(a^, *)!■ For the spin density 
Sy{x,t) the following result can be obtained 

Sy (x, t) 



exp(- 



(x — {hko/m — a)ty 
d2(l + 72f2) 



)- 
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(x — {hko/m + a)ty 
^2(1 + 72^2^ 
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According to this equation both pats of the initial 
wave packet moving along the x direction with differ- 
ent velocities are characterized by the opposite spin ori- 
entation (at the same time the average spin component 
Sy = J Sy{x,t)dr is equal to zero). 

Note that the components of wave function depend 
only on coordinate x, that leads to py = Py = 0, however 
the velocity Vy{t) ^ 0. Really, using the definition Eq.(4) 
it is not difficult to obtain 



Vx{t) 



hko 



Vy{t) = -asin(2fcoai)exp(-(-^)''). (17) 

As follows from these equations the average Vy velocity 
performs the oscillations in the transverse direction [zit- 
terbewegung or jittering) with the frequency 2koa and the 
damping time is determined by the parameter d/a. 



III. EVOLUTION OF TWO DIMENSIONAL 
PACKETS AT THE PRESENCE OF SPIN ORBIT 
COUPLING 

We consider now the evolution of two dimensional wave 
packet at the presence of spin orbit coupling. Let us 
consider the following form of the Gaussian packet at 
the initial moment t = 0: 



*(r,0) 



= m 



ipoxx/h) 



(18) 



where pox = ^ko is the average momentum and C = 
Xj^/i^d. Then, using a Green's fimction method we ar- 
rive after some algebra at the following equations for the 
components of spinor (in the momentum space) 



Ci(p,i) 



d , apt ip^t 
cos(— — I exp — - — -- 



2 j2 



p d 



(19a) 



^ . ,x d px + ipy . ,apt ipH 

C2(p,i) = i=r 2.sm(— )exp(-— — - 

^j-kTi p n 2mn 

P^d"^ _ fep^'^ , Pxkpd"^ 
2h' 2 ^ h 
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After that \l/i,2(r, i) can be obtained directly by 2D 
Fourier transform of Ci,2(r, t): 
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*i(r,i) = j exp 



i-i 



q^ht q^d^ kid'' 
2m 2 2 
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*2(r,i) 



{q^ kgd'^ -r'^ + 2ikQd^r cos <p) x 
X cos{aqt)qdq, 

id r cos y + sin -f — Ui[)d^ 
■v/tt -\/r2 - fcgd^ - 2ikod''rcosip 

q^M q^<e kl<P 



(20a) 



X j exp(— i 



2m 



X Ji(Y'(r2 - k^d"^ — 2ifcod2rcosiy?)x 

X ^va{a.qt)qdq, (20b) 

where J\ and /q are the Bessel and the modified Bessel 
functions of the first and the zeroth order, is asimutal 
angle in the xy plane. These expressions become simpler 
if the average momentum of a wave packet is equal to 
zero, i.e. pox = 0. In this case 



oc 

d f 

*i = — ^ / qJo{qr)cos{aqt)) 



.q^M q'd\ 

X exp(-i-- 7^)dq, 

2m 2 



(21a) 
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d y 



IX 



00 
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qJi{qr) sm.{aqt)x 



X exp(— f 



.q'^nt q^d^ 



2m 



)dq. 



(21b) 



As in the case of ID packet the shape of the full elec- 
tron density p{x, t) = \^i\' + 1*2 P at t > depends on 



the parameter 77 



In Fig. 2. wc show the elec- 



tron density p{x, t) for the case pox = at the time i = 5 
(in the units of ^) and r] = 2,7. As one can see the spin- 
orbit coupling qualitatively change the character of the 
wave packet evolution, so that during the time the initial 
Gaussian packet turns into two axially symmetric parts. 
As follows from our analytical and numerical calculation 
the outer part propagates with group velocity which is 
greater than a and the inner part moves with group ve- 
locity lower than a. If 77 <C 1, i.e. the packet is narrow 
enough, its evolution remained the standard broadening 
of the Gaussian packet of free particle. 

In Fig. 3(a) it is shown the packet evolution for the 
case pox = hko 7^ 0. It is clear that in this case the 
cylindrical symmetry is absent, and two maximums of 
the electron density spread along the x direction with 
not equal velocities. Each one of these two parts are spin 
polarized. Fig. 3(b) illustrates the distribution of the 
spin polarization Sy{x,y,t)for the initial state, polarized 
along z axis, Eq.(18). It is a smooth function which has 
different sign in the regions for two maximums of the 
electron density. 
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When Px ^ Q the motion of the wave packet center 
along X accompanied by the oscillation of the packet cen- 
ter in a perpendicular direction, or zitterbewegung. Below 
we consider the effect of damping of zitterbewegung oscil- 
lation for 2D packet which was not predicted in^. 

Using Eq.(22a) and Eq.(22b) we calculate the average 
value of the operator y — ifi-^ and obtain for t > the 
result 



2d2 



oo 

= -^exp[-(fcod)^] J sin2(^)x 





p^d\ ,2pkocP 
X exp( — -2-)-^i( — T — )dp. 



(22) 



In the case when wave packet is wide enough and the 
inequality a = dko ^ 1 takes place, one can obtain a 
simple asymptotic formula for y(t). To show this we rep- 
resent Eq. (22) as a sum of two terms 



velocity at the value of l/2fco- The last result coincides 
with*". Since the packet moves with constant velocity, the 
time oscillations of y{t) can be easily converted to the os- 
cillation of the wave packet center in real x, y space. 



IV. CYCLOTRON DYNAMICS OF 2D WAVE 
PACKET IN A PERPENDICULAR MAGNETIC 
FIELD 

In this section we examine the cyclotron dynamics 
of electron wave packet rotating in a magnetic field 
B(0, 0, B) which is perpendicular to the plane of 2D elec- 
tron gas. In this case the one-electron Hamiltonian in- 
cluding Rashba term reads 

(p + eA/cV ^ . . _ , 

H = h a[(jy[px + eAy/c) - (JxPy) + 
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(25) 



y{t) = —dexp{—a^) 



exp(— u^)/i {2au)du~ 



,2atu. , T, , ,„ , , 
cos( — ^) exp(— M )/i(2au)aM 



-dexp(— a^) 



where 



we denote \ = 



(cxp(a2) - 1) - Z 



Z = Re{ J exp(- 



(23) 



i^^j^)Ii{2au)du) . To evaluate Z we replace the modi- 
fied Bessel function Ii{2dkou) by it's asymptotic formula 
Ii{x) = , which valid for the case kod ^ 1. After 
that the integral with respect to u can be evaluated us- 
ing the stationary phase method that leads to the simple 
result 



d2 ' 



cos(2Qfcoi). 



Substituting this expression into Eq.(23) we finally have 



2ko 



1 — exp(^ 



) cos{2akot) 



(24) 



The last result demonstrates clearly that y{t) experi- 
ence the damped oscillations with the frequency 2ako 
decaying for the time ^. In the real 2D structures 
the frequency of the zitterbewegung have the order of 
10" - lO^^ggc-i foj. ^ iQ-5 _ io-6cni. The ampli- 
tude of the zitterbewegung is proportional to the electron 
wavelength in x. At Fig. 4 we plot the function y{t) de- 
termined by Eq. (22) which demonstrates in accordance 
with Eq.(24) the effect of zitterbewegung damping. When 
t ^ — the oscillations stop and the center of the wave 
packet is shifts in direction perpendicular to the group 



Here e is the electron charge, m is the effective mass, 
Px,y are the momentum operator components, a is the 
parameter of Rashba coupling, g is the Zeeman factor, 
and /^B is the Bohr magneton. Bellow we use the Landau 
gauge for the vector potential A = {—By, 0, 0). Then the 
eigenvalues and the eigcnfunctions of the Hamiltonian 
(25) indicating by quantum numbers n, kx, s = ±1 and 
corresponding to two branches of levels can be evaluated 
analytically (see, e.g.r^) 



2fc2 



huJrn ± En 



2na^h 



1/2 



(26) 



where E'^ = — gpsB is the zero Landau level, 
n = 1,2,3,..., LUc = is the cyclotron frequency, 

£b — \J is the magnetic length. The eigenspiniors 
are 



-iDn(j)n~l{y - Vc) 

e'''-"= / 0«-i(2/-yc) 



,J2ttA^ V "iDnipsiy - Vc) 







° ^/2^ \My~yc) 



(27) 



Here coefficients Z)„ are given by: Z?„ = 
i^o+V^+wnV^ ' A. - 1 + D^, 4>^{y - y,) are 
linear oscillator wave functions, j/c = ^%kx is the center 
of oscillator. It should be noted that for enough weak 
magnetic field the dependence of energy E^ on quantum 
number n (rt 3> 1) resembles the behavior of the function 
e_(p), Eq. (3). Namely, for small n the values of energy 
E^ are negative, decreasing with n, like for the hole 
states. 

Using the Eqs.(26),(27) we can obtain components of 
the matrix Green's function which permits us to find the 
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time evolution of the initial state. The usual definition 



Gij{r,r',t) = ^ dk^x 

oo 

(^C,fc.,i(r,i)C:fc.j(r,0). (28) 



n=0 



yields 



+00 



Gn(r,r',f) = i- J dA:,e^'=»(^-^') x 

— 00 

X X] fn+i{'t)^n{y - yc)4>n{y' - Vc), (29a) 



n=0 



where ip(x,y,7i) =iuj^ 2 4 2£b ■ 

The electron density obtained by numerical evalua- 
tions of the integrals Eqs.(32a) and (32b) is represented 
in Fig. 5 for relatively weak spin orbit coupling and 
strong magnetic field. The calculations was made for 
the material parameters of two dimensional GaAs het- 
erostructure: m = 0, 067mo, a = 3, 6 • 10^ cm-sec~^, g = 
-0,44, B =1T and k^x = Pox/h = 1,5 • 10^ cm-^ It is 
not difficult to verify that the series in Eqs. (32a) and 
(32b) converge very rapidly as n increases. So for our 
parameters it suffices to take Umax = 25 to calculate the 
components '(/'i(r,i) and '(/'2(r,t). At t > the initial 
wave packet (Fig. 5(a)) splits on two parts (Fig. 5(b)) 
which "rotate" with different incommensurable cyclotron 
frequencies. In accordance with Eqs. (26) these frequen- 
cies can be determine by the expressions 



+00 



G2i(r,r',i) = i- j rffc,e''=»(^-^')x 

— 00 

X ^9n+lit)(t)n+l{y - 2/c)^n(y/ - 2/c), (29b) 



E: 



n+l 



0Jc±\ El + 2{n+l)-;^T 



n=0 



TS^El + 2n^. 



(33) 



where the time-dependent coefficients fn{t) and gn{t) are 
given by 

fn{t) = e-''^="*(cos(5„i - i{l - ^) sinSnt), (30a) 



The effective n in this equation is connected with cy- 
clotron radius via relation RJt) = = \/2n^s- Be- 

lieving that <r = '^ITjl "C 1, i.e. in the case of a weak 



ff„(i) =e- ^ sin 



(306) 



and 5. 



2 , 2na^h^ 
i «5 • 



spin orbit coupling or strong magnetic field, one can ob- 
tain from (33) the approximate expression for the differ- 
ence between cyclotron frequencies 



(34) 



Eq 



Let the initial state coincides with the wave function 
of the coherent state in a magnetic field 



(31) 



Such choice of wave function 'I'(r, 0) is motivated by the 
following: as it is well known, at the absence of spin orbit 
coupling the dynamics of coherent states in a magnetic 
field looks like the dynamics of a classical particle. To an- 
alyze the time evolution in our case one needs to calculate 
the wave function at f > 0. Straightforward algebra with 
using Eqs. (29a), (29b), (31) leads to the final expressions 

' V2tt£b 2"n! J 



+ 00 

(r i) = ^ V / a 



x{-urHn+i{^-u) 



Fig. 5(b) demonstrates for the case <r <IC 1 the distribution 
of electron density at the moment when two parts are 
located at opposite points of the cyclotron orbit. The 
correspondent time can be determined from the relation 
(a;+ — io~)tQ = TT and hence for the GaAs structure we 
will have to = +^ _ = — tt4^ = 45^. 

After some cyclotron periods two split packets merge 
again which is demonstrated in Fig. 5(c). With time due 
to the effect of the broadening electron probability dis- 
tributes randomly around cyclotron orbit that is shown 
in Fig.5(d). 

In the opposite case of relatively strong spin orbit 

coupling or weak magnetic field when the inequality 
— 3> 1 holds true the difference between two 

cyclotron frequencies, as it follows from Eq.(33), equals 
to a;+ — u>~ = ^? . For the InGaAs structure with 
parameters m = 0,05mo, a = 3,6 • 10® cm • sec^^, g = 
— 10, B = IT and kox = Pox/^ = 1,5- 10® cm~^, we have 
? = 8 and the divergence time « 2, 3 • 

One can analyze the effects of the periodic splitting 
(|yev(2;,a,«) X a,nd reshaping of the wave packet in magnetic field as well 

as the process of distribution around cyclotron orbit by 
considering the time dependence of the cyclotron radius 
(^2^^) determined as R{t) = + {0)}^ ■ To do this we 



(32a) 
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represent the average value of coordinates xi = x, X2 = y 

as 

Xi = j ipi{r,t)xitl)i{r,t)dr+ 

+ J i>Ur,t)xi'i(>2{r,t)dr, (35) 

where -01 and ^2 are determined by Eqs. (32a) and (32b). 
The lengthy calculations (see the Appendix) eventually 
yield the explicit expression for the x{t) and y{t): 

x(t) = -%exp(-^2£^)x 

X {cos Wc* V] Sk{t)H2k+i (iJ \pox(-b/K)+ 
fe=o ^ ^ 

+ sm.ojct^Pk{t)H2k+i{i'\j^PoxiB/K)}, (36a) 



fe=0 



y{t) =q£% + — exp 



3 """P^ 3h' 



x{cosa;ctX]-^'^(*^-'^2fc+i(«Y -^Pox^b/K)- 
fe=o 

- sinwct ^ S'fe(t)iJ2fe+i(iy|^Pox^s/?i)}- (36b) 



fc=n 



As one can see, the dependence of x{t) and y(t) on 
the time are determined by both the factors cos Uct and 
sin Wei as well as by functions 

Sk{t) = i {-r7:)^[ik+2 cos4+iisin(5fe+2i- 
fc! 12 

-^fe+i cos 5k+2t sin 5k+i\ , (37a) 



Pkit) = i^^{—)''[cosSk+itcosSk+2t+ 



X sin 6k+it cos 6k+2t] , 



(37b) 



which describe the additional time dependence due 
to spin precession. Note that the frequencies Sk = 

Eq + '^^'^^^'^ are incommensurable. As a check on 

this formalism, it is not difficult to show that in the ab- 
sence of Rashba coupling {a = 0) as it follows from Eqs. 
(36a) and (36b) 



x{t) = poxi%/hsmoJct, 

y{t) =POxi%{^ - COSWct), 



that correspond to the classical motion of charged parti- 
cle in the magnetic field with a constant radius. 

The time dependence of the cyclotron radius R{t) in 
the system with Rashba is presented at Fig. 6. It is 
clear that the oscillations of R{t) are connected with the 
effects of periodic splitting and reshaping of wave pack- 
ets. The radius has the minimal values at the moments 
when two parts of the packet are located at the oppo- 
site point of cyclotron orbit. This situation is shown at 
FIG 5.(b). The first minimum labeled by the letter b at 
FIG .6 One can see that the time of the first minimum 
approximately coincides to our estimation made above: 
to ~ 4:bTc. The radius is maximal at the moments of 
the packet reshaping that is shown at Fig. 5(c) (two of 
these points labeled by the letters a and c. Due to the 
effects of incommensurability of the cyclotron frequencies 
and the packet broadening the amplitude of the oscilla- 
tions decrees with the time After that the electron density 
distributes around cyclotron orbit, the amplitude of the 
oscillation ceases and the electron density distribution 
acquire the no regular character (Fig. 5(d)). 

We evaluate also the distribution of the electron den- 
sity for the structure with relatively strong spin orbit 
coupling. For such systems instead of the repeated pro- 
cess of the splitting and restoring of the wave packet 
discussed above the transition to the irregular distribu- 
tion along the cyclotron orbit can be realized for the 
time of the order of one cyclotron period. This conclu- 
sion is confirmede by the simple estimation made for the 
InGaAs/GaAs structure discussed above. 



V. CONCLUSIONS 

We have analyzed the evolution of ID and 2D wave 
packets in 2D electron gas with linear Rashba spin orbit 
coupling. We showed that the electron packet dynamics 
differs drastically from usual quantum dynamics of elec- 
trons with parabolic energy spectrum. Depending on the 
initial spin polarization packet splits in two parts which 
propagate with different velocities and have different spin 
orientation. At the time when two parts of wave packet 
overlap, the packet center performs oscillations in much 
the same way as for a relativistic particle. The direction 
of these oscillations is perpendicular to the packet group 
velocity. When the distance between split parts exceeds 
the initial width of the packet these oscillations stop. 

In the 2D semiconductor structures placed in a perpen- 
dicular magnetic field the spin orbit coupling changes the 
cyclotron dynamics of charged particles. As at the ab- 
sence of magnetic field the initial packet splits in two 
parts, which rotate in a perpendicular magnetic field 
with different incommensurable cyclotron frequencies. 
As a result, after some cyclotron periods these parts join 
again. The corresponding time to essentially depends 
upon the ratio of the energy of spin orbit coupling and 
the distance between Landau levels, Eq.(26): <; = '^^rj^i ■ 



(38) Thus, for the systems with weak and relatively strong 
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spin orbit coupling e.g. GaAs and InGaAs hetrodtruc- 
tures, the time equals to 45Tc and 2, STj,, respectively. 

The splitting and zitterbewegung of the wave packets in 
nanostructures with spin-orbit coupling can be observed 
experimentally in low dimensional structures. In particu- 
lar, these effects should determine the electron dynamics 
and high-frequency characteristics of the field effect tran- 
sistor by Datta and Das^'^, and other spintronic devices. 
Simple estimations show that during the time of the wave 
packet propagation through the ballistic transistor chan- 
nel where the distance between emitter and collector is 
of the order of Ifim, the distance between two split parts 
of the wave packet becomes comparable with its initial 
size. In this situation the high-frequency characteristics 
of the field effect transistor should be substantially af- 
fected by the spin-orbit coupling. Moreover, the atypical 
semiclassical dynamics of a spin-orbit system placed in 
a magnetic field will influence the shape of the cyclotron 
resonance line in 2D systems with spin orbit coupling. 
An important feature of these experiments is that the 
electron transport is in the ballistic regime and thus the 
momentum relaxation time t should be considered much 
more greater compared with the typical splitting time. 
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VI. APPENDIX 

This appendix provides some of details involved in ob- 
taining the average value of the position operator given 
by Eqs.(36a), (36b). According to Eq.(35) 



functions 



Consider the calculation of the first term yi{t). Using 
the initial wave function, Eq.(31), we obtain 

dr'dr" 



X exp 



21% 



X J Gn{r,r',t)yGl,{r,r",t)dr. (A.2) 

Denote the last integral in this equation as 

Mfi- I Gn{r,r\t)yGl,ir,r",t)dr. (A3) 

Then substituting Eq.(29a) into Eq.(A.3) and using 
the well known formula for a linear harmonic oscillator 



+ 00 



v4>n{y - yc)4>k{y ~ yc)dy 



^ \^Vn5k,n-l + VnTlSk-n+ij + yc5nM, (A. 4) 



we will have 



+ 00 

-hi' 



^ik^{x-x")fi{y' ,y" ,t.,ya) 



{A.5) 



Here 



X ( E V^fn+i{t).f:{t)My' - 2/c)0„-i(y" - yc)+ 

+^v^/„+i(t)/,:+2(t)x 

n=0 

x0n(y' - yc)(/>«+i(y" - yc)^ + 
+yc J2 \fn+i{t)\^My' - yc)My - yc), (A.6) 

n=0 

where the coefficients /„(t) are given by Eq.(30a). We 
calculate yi{t) by substituting Eqs.(A.5), (A.6) into 
(A.2). The resulting integrals can be evaluated by us- 
ing Gaussian transformationi^ 

^ J e-(--^)'iJ„(2/)dy=(2x)", 

Hniix) 



— OO 
+ 00 



(2i)" 



(A.7) 



{A.l) Finally we obtain 



yiit) = — exp( 



6 """^ 

^Mt)H2k+l{tVV^P0cciB/h), (A., 



k=0 



where 



+/fe+iW/fc+2W-2iA+i(on- 



(A.9) 



Performing the same kind of calculation we have for y2{t)'- 
y,(i) = -exp(-^^)x 
x^7fc(t)i/2fe+i(«y273pox^B/?i), (A.IO) 

fc=0 
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where 



7fcW 




i9*k+iit)9k+2{t)- 



-9k+i{t)gl+2{i) - 2|gfc+i 



(A.ll) 



and the coefficients fk{t) and gk{t) in Eqs.(A.9), (A. if) 
are determined by Eqs.(30a), (30b). The preceding ex- 
pressions immediately lead to the average value y{t) given 
in Eq. (36b). The evaluation of x{t) can be obtained 
by following the procedure similar to that which led to 
Eq.(36a). 




FIG. 1: The electron probability density p[x,t) = + 
|5'2p(a), and spin density Sz(b). The dashed, thick and thin 
lines correspond to different moments of the time namely ti — 
0, t2 = 1, 5 and — 7(in the units to = 7^^)- 




FIG. 2: (Color on line). The electron probability density 
p{x,t) = l^'ip + 1^2p for the initial state Gaussian packet 
Eq.(18) with pox = at the time t = 5 (in the units d/a). 
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FIG. 3: (Color on line). Electron density p{x,t) = I'Jfif + 
|>I'2p(a) and spin density Sz{x,y,t) (b) for kgd = 1 at the 
moment f = 5 (in units d/a). 




FIG. 4: The average coordinate of y{t) versus time for the 
packet with kod — 25. 



FIG. 5: (Color on line). Evolution of coherent wave packet 
Eq.(31) in a perpendicular magnetic field: (a) the initial 
electron density Eq.(31), (b) two split packets at time to ~ 
452p,(c) restored packets at time 2to ~ 90^, and (d) ran- 
domized electron density for large time. 




'/Tc 



FIG. 6: Cyclotron radius plotted versus the time (for the same 
parameters as in Fig. 5) . The distance between maximums and 
minimums of R{t) marked by arrows is approximately equal to 
SOTc Time is measured in units of cyclotron period 7c = 
The points a, b, c, d correspond to the same moments of time 
as in Fig. 5(a), (b), (c), (d), respectively. 



